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We consider three-dimensional reshaping of thin nemato-elastic sheets containing half-charged defects upon 
nematic-isotropic transition. Gaussian curvature, that can be evaluated analytically when the nematic texture is 
known, differs from zero in the entire domain and has a dipole or hexapole singularity, respectively, at defects of 
positive or negative sign. The latter kind of defects appears in not simply connected domains. Three-dimensional 
shapes dependent on boundary anchoring are obtained with the help of finite element computations. 


Liquid crystal elastomers and glasses (T), made of cross- 
linked polymeric chains with embedded mesogenic structures, 
combine orientational properties of liquid crystals with shear 
strength of solids, and were envisaged by de Gennes as pro¬ 
totype artificial muscles 0. Their specific feature is a strong 
coupling between the director orientation and mechanical de¬ 
formations, which can be controlled by the various physical 
and chemical agents. When the material undergoes a phase 
transition from the isotropic to nematic state, it strongly elon¬ 
gates along the director and, accordingly, shrinks in the nor¬ 
mal direction to preserve its volume; the opposite effect takes 
place as result of the reverse transition. Stresses arising due 
to these intrinsic deformations were investigated for flat sheets 
where they were shown to lead to phase separation of isotropic 
and nematic domains 0 or formation of persistent defects 
that are not necessitated by topology l4l . 

Internal stresses can be relaxed when three-dimensional 
(3D) deformations are allowed. The reshaping effect causes 
bending of flat thin sheets into curved shells. Aharoni et al (5) 
studied the various nematic textures that can be impressed in 
the material to induce a desired two-dimensional (2D) metric 
that would determine a 3D shape upon transition from nematic 
to isotropic state (NIT). This opens the way to construct sur¬ 
faces with nonzero Gaussian curvature out of flat sheets with 
no stretching m, or, more generally, to change the Gaussian 
curvature of shells. While the Gaussian curvature is an intrin¬ 
sic property that can be computed based a surface metric only 
q, computing an actual 3D shape is a far more challenging 
task, which is still more complicated by constraints on ne¬ 
matic tissues imposed by boundary conditions or topology of 
closed manifolds that necessitate emergence of defects. A 3D 
shape can be computed analytically only for a symmetric tex¬ 
ture, and the only available example is deformation of a circle 
containing a defect with unit circulation into a cone or an an¬ 
ticone mm. This case is exceptional in two respects. First, 
the natural nematic texture induced within a circle by either 
tangential or normal boundary anchoring contains two half- 
charged defects which have a lower energy than a single defect 
of unit charge. Second, the Gaussian curvature induced by a 
unit charge defect vanishes everywhere except the vicinity of 
the defect itself s, where it can be regularized either by local 
stretching m or by the depletion of nematic order 0. Away 
from the defect, where the Gaussian curvature vanishes, fold¬ 
ing turns out to be necessary 0 to accommodate an extended 
perimeter of the circle into a sphere with a shrinking radius, in 
a way similar to crumpling of isotropic non-extensible sheets 
0. 


In a more difficult and realistic case of half-charged defects, 
even planar textures can be computed analytically only in sim¬ 
plest configurations, and the lack of axial symmetry leaves lit¬ 
tle hope for analytical theory of 3D shapes. The Gaussian cur¬ 
vature differs then from zero also outside defect cores, which 
makes smooth 3D shapes possible. In this communication, 
we compute numerically the shapes generated upon NIT from 
naturally formed nematic textures. We first summarize our 
basic methods: analytical, starting with metric transforma¬ 
tion induced by NIT and extending, in the case when the ne¬ 
matic texture is known, to computation of Gaussian curvature, 
and numerical, indispensable for visualizing actual 3D shapes. 
We consider next representative particular cases: a vicinity of 
isolated defects, where the induced structure can be well re¬ 
solved; an elliptic sheet containing two half-charged defects; 
and a sheet with two holes where negatively charged defects 
naturally appear. 

We consider a thin sheet with preferred parallel orientation 
of the director on both outer surfaces. The material is de¬ 
scribed then by the two-dimensional 2D nematic order param¬ 
eter, presented as a traceless symmetric tensor dependent on 
2D coordinates: 


Q a p = S (n a np - \Sfj , ( 1 ) 

where n a are vector components of the director, 5% is the Kro- 
necker delta, and S is the scalar order parameter equal to unity 
in a perfectly aligned nematic state. Under the assumption 
of negligible elastic deformations, an infinitesimal interval is 
transformed as a result of a change Q a p of the nematic state 
as 


dr=TVA T a p = TV- 1 / 2 (^ + aQ a fj , (2) 

where a is a metric factor measuring the relative extension and 
contraction, and TV is a normalization factor; summation over 
paired indices is presumed. The original metric tensor g° is 
transformed thereby to g = Tg°T; the transformation can be 
made area-preserving 0 by choosing N to ensure det(T) = 1. 
In the case of transition from nematic (S = 1) to isotropic 
(S = 0) state (NIT), the explicit metric expression is 

9a p = (1 - a 2 )” 1 [(1 + a 2 )g% - 2 aQ a0 ). (3) 

If the coordinate axes x a , are oriented, respectively, along 

and normal to the director, Eq. © simplifies to 


d^ 1 = ^dx 1 , dp=£dx 2 , 


l-a 


( 4 ) 
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and the transformed metric is diagonal with the elements 

r 2 , e. 

The transformed metric can be used to compute intrinsic 
properties of bent sheets or shells but not their extrinsic prop¬ 
erties defining their position in space. Among the former, the 
most important one is Gaussian curvature equal to the only 
non-zero component K = R 1212 of the Riemann curvature 
tensor 0 . It is defined by standard formulas dependent on 
the first and second derivatives of the metric tensor. 

Computing actual 3D shapes which, unlike Gaussian cur¬ 
vature, are not determined by intrinsic characteristics of a sur¬ 
face alone, is a much more difficult task that can be only ac¬ 
complished numerically. We carried out computations by tri¬ 
angulating the surface and defining the nematic order param¬ 
eter at the grid nodes. The nematic director on triangular tiles 
was defined as the average of the three vertices. Normal vec¬ 
tors were computed for each tile, and the normal vectors 
at nodes mj were defined as the average of the normals at the 
surrounding tiles. If the nematic field was known analytically, 
as in some following examples, ??, the nematic director was 
defined at the centres of triangular tiles. 

The triangles were addressed at random following a Monte 
Carlo routine, and NIT was imitated by changing the length of 
their sides (edges) according to Eq. ([2]) and further rescaling to 
eliminate numerical errors in maintaining area conservation. 
The procedure continued, typically, for 20 Monte Carlo cycles 
to minimize the energy functional defined as 

E = ~ l°i) 2 + K ' m *i)- 

edges nodes j 

Here, the first term accounts for deviation of the length li of 
an ith edge from the value defined by Eq. ([2]). In the sec¬ 
ond term, m^ denotes the normal to the a j th tile of those 
surrounding an it h node. The sum over all neighbors mea¬ 
sures the deviation from their average orientation mwhich 
accounts for the curvature of the bent shell. Since bending 
rigidity is much weaker that in-shell shear modulus, the ratio 
n of bending rigidity to in-shell elasticity should be small. In 
order to speed up numerical algorithm at early stages, it was 
gradually decreased from 10 -4 to 5 x 10 -6 , which allowed 
us to arrive as close as possible to the required degree of local 
extension or shrinkage. 

When the nematic texture was not known analytically, we 
computed it in one-constant approximation by minimizing the 
nematic energy f(Vn) 2 dE over a studied domain with appro¬ 
priate anchoring conditions at boundaries. The derivatives are 
approximated by differences between the directors at neigh¬ 
boring nodes. Although the director is more naturally defined 
on tiles, redefining it on nodes improves numerical accuracy 
due to a larger number of neighbors (six rather than three). 
We resolved the defect core structure only in the case of iso¬ 
lated defects leaving otherwise the scalar order parameter S 
equal to unity, but triangulation density was increased in the 
vicinity of defects. An amended procedure is applicable to 
computation of textures on curved surfaces. In this case, co¬ 
variant derivatives are computed by virtually rotating neigh¬ 
boring tiles into a common plane before comparing the direc¬ 
tors, and an additional energy penalty is introduced according 


to the surface curvature along the director. 

Strong reshaping effects in the vicinity of defects are of 
special interest, and can be studied by considering isolated 
defects, far removed from other defects and boundaries. The 
texture in the vicinity of positive and negative half-charged de¬ 
fects has different symmetry (Fig.|T]), even though the nematic 
order parameter is expressed in a similar way. In a vicinity of 
defects of the charge ±1/2 (but outside their cores) the tensor 
order parameter is 


Q = 


1 

2 


cos (j) ±sin0 \ 
± sin 0 — cos 0 ) 5 


( 6 ) 


where 0 is the polar angle. Starting from the Cartesian metric 
g°, we obtain using Eq. ([3} the transformed metric 

_ 1 f 1 — 2a cos 0 ± a 2 =f2 a sin 0 

^ 1 — a 2 \ ±2a sin 0 1 ± 2 a cos 0 ± a 2 

( 7 ) 

The Gaussian curvature K± of a bent shell formed upon NIT 
is derived from this metric: 

acos0 3a cos 30 

= ~2n -2± = —oTl-2V (°) 

r z { 1 — a z ) r z { 1 — a z ) 

These simple expressions are in accordance with the symme¬ 
try of the bent shell around the defect. The curvature diverges 
at the defect location (r —0) leading to a dipole or hexapole 



FIG. 1: The vicinity of isolated positive (left) and negative (right) 
half-charged defects. Above: the director field; middle row: the reg¬ 
ularized Gaussian curvature; below: 3D shape. 
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singularity, respectively, for defects of positive or negative 
sign. For the positive charge, the curvature is positive near 
the defect’s “comet tail” at = 0, while a saddle geometry 
with negative Gaussian curvature is indicated on the opposite 
side. For the negative charge, there is a threefold symmetry, 
matching the symmetry of the nematic field (Fig. [TJ. 

The divergence can be regularized by taking into account 
vanishing S at the defect core. For defects of either sign, the 
radial dependence S(r) is identical to that of the amplitude of 
a complex scalar field near a defect of unit charge ED, and 
satisfies the equation 

S"(r) + r _1 S"(r) + S (l - r~ 2 - S 2 ) = 0, (9) 

where the radial coordinate is scaled by the healing length. 
The regularized Gaussian curvature within a circle encom¬ 
passing the defect core is shown in Fig. |T] 

In only two known examples, nematic texture containing 
half-charged defects can be obtained analytically in the one- 
constant approximation by means of conformal transforma¬ 
tion Col. The simplest case is an ellipse where the director 
is aligned along coordinate lines u = const (if it is anchored 
normal to the boundary) or v = const (if the anchoring is tan¬ 
gential) of the elliptic coordinate system, which are identical 
to the isolines of the real and imaginary part of the analytic 
complex function w = sin 2 (Fig. [2]). Upon NIT, the metric 
of the resulting curved sheet is readily obtained by modifying 
the elliptic metric according to Eq. 0: 

g = gdiag (/ _2 ,/ 2 ) , g = ^(cosh2u — cos2 u), (10) 

where the half-focal distance is normalized to unity. This 
expression is equally applicable to textures with the director 
along the coordinate lines u = const or v = const; one needs 
only to interchange the metric factors corresponding to the 



FIG. 2: Above: director orientation for normal and tangential bound¬ 
ary anchoring. Below: Gaussian curvature. 



FIG. 3: Deformation of an ellipse. Above: 3D shapes for normal an¬ 
choring and tangential anchoring (a = 0.1). Below: the dependence 
of the mean square deviation of the surface from the original position 
as a function of the metric factor a. 


two variables. In both cases, the Gaussian curvature is com¬ 
puted as 


K = ±- 
9 

= ± 


<9/1 dg\ <9/1 dg 
du \g du) + dv \gdvj 
16a(cosh 2v cos 2 u — 1) 

(1 — a 2 )(cosh2v — cos 2 u) 2 ’ 


(ID 


This expression has a dipolar singularity at both foci (see 
Fig. |5J. Near the foci at x = ±1, y = 0, one can set 
u « ±^/2r sin(</>/2), v ~ :±i/2r cos(0/2) and recover the 
first formula in ([8]) in the leading order of expansion in the 
distance r <C 1 from either focus. The difference is in the 
position of “comet tails” of both defects, which are directed 
toward the centre of the ellipse when the anchoring is tangen¬ 
tial, and outwards when the anchoring is normal. The pos¬ 
itive and negative curvature regions are situated accordingly, 
so that the Gaussian curvature plot in Fig. [^corresponds to the 
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normal anchoring and has to be flipped over when the anchor¬ 
ing is tangential. The shapes for both cases shown in Fig. [3] 
are quite dissimilar. The lower panel in Fig. [3] shows the non¬ 
linear dependence of the mean square deviation of the surface 
from the original position as a function of the metric factor a. 

Negatively charged defects appear when domains are not 
simply connected. Inserting a single hole with the same an¬ 
choring as on the outer boundary leads to a lowest-energy tex¬ 
ture with no defects. When a second hole is added, two de¬ 
fects with the charge —1/2 must appear. The position of de¬ 
fects depends on the size and position of holes. Some typical 
examples of elliptic domains with symmetrically placed two 
holes, as they deform following NIT, are shown in Fig. [4] In 
all cases, the anchoring at all boundaries is normal. 

In the upper left picture with relatively small holes, the de¬ 
fects are situated symmetrically on the line connecting the 
centres of both holes. In the upper right picture, where the 
holes are large, the defects are placed on the symmetry axis 


separating the holes. The arrangement of defects in the two 
lower pictures with holes of intermediate size is less intuitive. 
In the left picture, the defects are situated asymmetrically near 
one hole, while in the right picture, they are placed diagonally 
near alternative holes keeping a larger mutual separation. The 
latter arrangement has a lower energy, so that its counterpart 
is metastable. The corrugated shapes obtained following NIT, 
are, naturally, different in all cases. 

Reshaping of nemato-elastic sheets or shells opens ways 
of creating a variety of forms that can be manipulated by 
boundary anchoring, positioning of defects, and topological 
changes. Besides static reshaping, the forms can be actuated 
dynamically, thereby creating crawling and swimming micro 
robots that will be the subject of a further study. 
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FIG. 4: Deformation of domains with two holes (a = 0.1). Small 
circles mark positions of defects with the charge —1/2 
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